Acoustic attenuation in a type-II superconductor at high magnetic fields 
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We have calculated the longitudinal acoustic attenuation in a type II superconductor in high 
magnetic fields within a mean-field BCS theory. We predict two new features in the corresponding 
attenuation signal as compared to that of the Meissner state. Our analytical calculations predict 
the existence of oscillations in the attenuation as the external magnetic field is varied- this effect is 
associated with the Landau level structure of the electron states and is analogous to the well-known 
de Haas van Alphen oscillations in the mixed state. The attenuation directly probes the quasiparticle 
energies; the presence of gapless points in the quasiparticle spectrum, which is characteristic of 
type-II superconductors at high magnetic fields, shows up in the frequency lo and temperature T 
CI dependence of the attenuation in the limit of low uj and and low T respectively. At low T there is 

no analogue to the discontinuity in the attenuation observed in the Meissner state when Tiuj — 2 A, 
where A is the quasiparticle energy gap. This result opens up the possiblity of experimentally 
Q ■ determining the existence and nature of the gapless points in the quasiparticle spectrum of a type-II 

O ' superconductor in high magnetic fields. 

CLf PACS numbers: 74.25.Ld, 74.60, 
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I. INTRODUCTION 



The problem of understanding type-II superconductivity in high magnetic fields has been a lively topic for some 
years now. When fcgT < Tilu c , where uj c = eH/mc is the cyclotron frequency, the Landau level quantization of 
the electron levels become important. This gives rise to a number of interesting effects.El One such effect is the 
existence of magnetic oscillations in the free energy of the_superconductor in the mixed state. These de Haas-van 
Alphen (dHvA) oscillations was observed over 20 years agca and have since been observed in a number of materials. 
One important theoretical consequence of the Landau level ^quantization is the presence of gapless points in the 
quasiparticle spectrum. Theoretical investigations have showrtra that the quasiparticle spectrum of the mixed state 
in high magnetic fields is characterized by a set of gapless points in the magnetic Brillouin zone (MBZ). These 
gapless point leads to an algebraic temperature dependence of thft-thermodynamic functions and an algebraic voltage 
' dependence in the tunneling conductance.!] It has been suggested^ that the existence of the dHvA oscillations in the 
mixed state is a consequence of these gapless points. It turns out however, that it is rather difficult to develop a simple, 
y— I ' yet consistent theory for the dHvA-oscillations since both the oscillatory behaviour of the ground state energy and 
the gapless nature of the quasiparticle energies need to be taken into account Jl The interpretation of the experimental 
results for the dHvA-oscillations is consequently somewhat unclear and does not, in our opinion, give a completely 
unambiguous signature for the presence of gapless modes. 

In this paper we consider the attenuation of longitudinal acoustic waves in the mixed state where the order parameter 
is assumed to form a vortex lattice. Since the absorption of the phonons is due to quasiparticle excitations, the 
experiment directly probes the quasiparticle density of states. For clean materials, theoretical results have proved 
rather difficult to obtain since an expansion in powers of the order parameter does not converged and the semiclassical 
approximation leads to unphysical results. □ Since we consider the superconductor at high magnetic fields, it is crucial 
O ■ to take into account the Landau level quantization of the electronic levels. It turns out that by including this effect, 
one avoids the difficulties encountered in the semiclassical approximation as predicted by Scharnberg.ta We calculate 
analytically the attenuation in two limits: ksT <C huj and fc^T 3> huj where u is the frequency of the sound wave 
and T is the temperature. We predict that the attenuation will be an oscillatory function of the external magnetic 
field due to the Landau level structure, in analogy with the usual dHvA oscillations. The frequency and temperature 
dependence of the attenuation is determined by the existence and the nature of the gapless points in the quasiparticle 
spectrum for fcgT <C hu and fcgT hu> respectively. This in principle gives an experimental tool for probing 
the nature of the quasiparticle energies in the mixed state. Our analytical theory is supported by essentially exact 
numerical calculations for the attenuation. 

As the analysis in this article will show, the existence of the gapless points in the quasiparticle spectrum has 
significant consequences on the longitudinal sound attenuation. We would therefore expect that this should also 
hold for systems with intrinsic gapless points even for no external field such as high-temperature superconductors 
(high-T c 's), which seems to have a cZ-wave gap symmetry. Hence a similar analysis of the acoustic attenuation in 
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these high-T c 's as the one presented in this paper for type-II superconductors in high magnetic fields would be very 
interesting, and we will return to this in a future publication. 

II. GENERAL THEORY 

We consider a weak-coupling superconductor in three dimensions (3D) described within mean field theory by the 
following Hamiltonian: 

^ = E/ <fr^(r) ( (P ~| A) - M ) Vv(r) + J dr[A(r)^(r)Vl(r)+cc] (1) 

where the order parameter is defined as A(r) = (/(^(r^j^r)) , g is the coupling strength and /i is the chemical 
potential. Our theory describes longitudinal ultrasonic waves in clean samples for high magnetic fields. We have 
for simplicity confined our theory to the case when q || H where q is the wavevector of the sound wave and H is 
the magnetic field. This simplifies our calculations since the coupling between the collective modes associated. with 
fluctuations of the order parameter and the longitudinal phonons can be neglected to a good approximation.Bii3 The 
attenuation of sound waves is given by the imaginary part_qf the retarded phonon self energy. To lowest order in the 
electron-phonon coupling we get for the attenuation a(w)tJ 

a(w) cx -Lolm {L> fl (q, u)} (2) 

where D R (q, oj) is the retarded density-density correlation function 

iD R (x, x', t-t') = ([n(x, t), n(x', t')])6(t - t') (3) 

and n(x, t) = n(x, t) — (n(x, £)) is the operator describing density fluctuations. This function is evaluated as the 
analytical continuation of the thermal correlation function. We will in this paper for simplicity ignore the electronic 
Zeemap effect. The effect of a finite spin splitting on the magnetic oscillations is well understood both in the normal 
stateE3 and in the mixed state.B Hence ignoring the spin splitting we obtain 

D R (q,u) = V / d 2 r [ d 3 r'e^' [G(r, r', u v )G{v' , r, w v - w 7 ) 

Vcell ~ J cell J 

-Ft(r, r', Uv )F(r>, r, Uv - u> 7 )] \^ +lS (4) 

where G(r, r', lo v ) and F(r' , r, oj v ) are the one-particle Green's functions for the superconductor. The symbol J u d 2 r 
implies integration over one vortex lattice cell in the x—y plane whereas J d 3 r' means integration over the whole crystal. 
By treating the electron-phonon matrix element as an overall factor in our formalism (i.e. its frequency dependence 
u> is included as a prefactor in Eq.(|2j)) we .have made use of the fact that the screening for longitudinal modes is 
essentially the same as in the normal phase.Ej Since the order parameter is assumed to form a vortex lattice, we can 
use a set of single particle states (r) characterized by a Landau level index n, a wavevector k = ik x , k y ) in the 

MBZ, and a wavevector k z along the z-direction. These basis states were introduced by Norman et alB In this basis 
there is no mixing between different k's since the sound wave travels along the z-direction. Expanding the Green's 
functions in this basis and using the k-space symmetry of the problem due to the presence of the vortex lattice,0 we 
obtain the following expression for the ultrasonic attenuation in the mixed state of a type-II superconductor: 

a(q, ^^j^-YsY.! dk *i [/(^ - f( E 'W(E' -E- u)(\U' A U\ 2 ~ 

X y nn' k 

—V' A V* U'* At/) — 5{E - E' - uj){\V* A V\ 2 - U'* A UV' A V*)] 

+[1 - f(E') - f{E)]5(E' + E- w)<\U' A V\ 2 + V' A UU'* A V*)} . (5) 

Here U = U£ k ^(r), V = V£ k ^(r), U' = U£ kr+q (r) and V = Vkk,+ q ( r ) are tne Bogoliubov functions. The quasiparticle 
energies are given by E = E^ (k) and E' = E k ^ +q (k); f(E) = l/(exp(E/kBT) + 1) and the A-product means 
integration over the ccy-plane. The quasiparticle energies, and associated wavefunctions can be found by solving the 
corresponding Bogoliubov-de Gennes (BdG) equations.^ We have set up a program that solves the BdG equations 
self-consistently in 3D such that we get the Bogoliubov functions as a function of H. The method for solving these 
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equations numerically in 2D is described in detail elsewhere .tfflThe extension to 3D is straightforward although it 
is computationally significantly more demanding. We choose the normal state dispersion law along the z-direction 
to be either the plane wave form e(k z ) = k 2 /2m, or more suitable for layered structures, the tight-binding form 
e(k z ) — tcos(k z a z ) where a z is the distance between the planes. 

In order to develop an analytical theory for the attenuation we need to make some approximations. Near the upper 
critical field H C 2 and for quasiparticle levels close to the Fermi energy we can, as a first approximation, igware the 
off-diagonal pairing (i.e. the so-called diagonal approximation). The quasiparticle energies are then given byE-3 

££(k) = V£n(M 2 + |A(k)| 2 (6) 
and the corresponding Bogoliubov functions are given by 




Here £„(fc z ) = ntuv c + e(k z ) — (i. It should be noted that this approximation is only valid for energy levels close to the 
Fermi energy. Further away from the Fermi level there are degeneracies between electron and hole states belonging 
to different Landau levels, hence our approximation will eventually breake down. However, for low frequencies 
and temperatures one can show that only levels for which Eq.(^|) and Eq.(^) hold will contribute to the damping 
described by Eq.(|5|). A closer examination shows the requirement for the diagonal approximation to hold to^ba 
max(fcsT, wLpJkj c /4. This follows because mixing with hole levels become important for levels with E > Huj c /4:Bt3 
Dukan et amii3 argued that the off-diagonal pairing does not change the qualitative behavior of the superconductor 
in a high magnetic field for fields not too far below H C 2, and that the quasiparticle spectrum remains essentially the 
same when the off-diagonal terms are included. From the diagonal approximation, it directly follows that there are 
gapless points in the MBZ (i.e A(k) = 0). Even when the diagonal approximation breaks down, there will still be 
points in the MBZ where the gap vanishes. It seems the role of the off-diagonal terms is to shift the value of tha 
Fermi momentum k z j where the gapless behaviour occurs away from its diagonal approximation value e(k z t) = O.u 
Eventually, true gapped behaviour sets in when the. superconducting order is strong enough to increase the energies 
of the quasibound states in the vortices above huj c B 

In Fig. 1 we have plotted the lowest quasiparticle energy, calculated numerically, along the T—M direction in k-space 
for two different values of the order parameter A(r) at a low temperature. The two points T = (n/4a x , ir/2a y ) and 
M = (3tt /4q x ,Tr/2a y ) are two of the corners in an irreducible triangle reflecting the symmetry of the BdG-equations 
in k-space.u Here a x — 1(\^3tt/2) 1 / 2 as we have chosen a triangular symmetry for the vortex lattice, a y = l 2 /a x 
and I 2 = Tic/eH is the magnetic length. We have chosen the size of the system such that ir/(2a x Ak x ) = 50 where 
Ak x = 2tt / Lx and L x is the extend of the system in the x-direction. The dashed lines give the diagonal approximation 
to the energies and the solid lines the exact numerical result. The two highest-lying curves are calculated with 
(A(k)) ~ 0.3?iw c and the two lowest lying curves are with (A(k)) ~ 0.05hui c . Here (A(k)) means the k-space average 
of the diagonal matrix element in the BdG equations. There are 10 Landau levels within the pairing width and 
the dispersion law along the z-direction is e(k z ) = t cos(k z a z ). The value of k z is chosen such that e(k z ) = 0. As 
can be seen, the diagonal approximation predicts two gapless points along this k-space direction, both with a linear 
dispersion law. For a small pairing parameter there is good agreement between the diagonal approximation and the 
full calculation. For a larger pairing parameter, i.e. (A(k)) ~ 0.3?Lu; c , the approximation is less precise- indeed the 
gapless points predicted by the diagonal approximation disappear in the full self-consistent calculation. This is due 
to the off-diagonal pairing which become increasingly important as the pairing interaction increases deeper into the 
mixed state. The gapless points are now at a different value of k z . We expect our theory to be valid reasonably close 
to the transition line such that we can ignore the off-diagonal pairing. As will be shown later, the major contribution 
actually comes from the gapless points where the diagonal approximation is most valid. We will now calculate the 
attenuation in two limits using this approximation. 



III. LOW FREQUENCY 



We will first treat the case uj <C fcsT. In this limit we can focus on the first two delta functions in Eq.(||) which 
describe the scattering of a quasiparticle. We will outline the calculation for the dispersion law e(k z ) = k 2 /2m. 
Making the approximation £'/E' ~ £/E which is valid for hui <C fcsT the first two terms of Eq.(||) become: 

a(q,w) = "^^^ / dk ^-K(n,k))d E f(El(k)) 1 -^^ (8) 
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Here k*(n, k) is the solution to the equation Ejt +g (k) — E% (k) = to. From the Poisson formula we have the identity 

^/(^(k))|^= £ fd^-d E f(Et t (k))M^ (9) 



E 



where = 35^c + e (^z) — A* = I 2 - — nf)huj c + e(&*). Here u/ = fj,/huj c — 0.5 is the Landau level index at the 

chemical potential. We have extended the Landau level sum to go from —00 as the low levels do not contribute 
anyway. Making the variable substitution £ x (k*) = zhuj c we end up with 

2 ®® /* OO I y- I 

a(q,c) = --f^E E e 2 ™^-<^/^) / dz e ^ fe/ (£ 2 ( k ))JS£L (io) 

where £ 2 (k) = (£ 2 + A(k)) 1 / 2 . We have approximated k* as the corresponding normal state solution (i.e. it is 
independent of n) . This approximation which is only necessary when we calculate the oscillations of the attenuation 
(i.e. the terms in Eq.(|To|) with j 7^ 0), should be good close to the gapless points and is equivalent to ignoring any 
phaseshift in the oscillations due to the superconducting order. The O'th harmonic is given by the j = term. We 
get 

2oj 2 m \ - 1 
a{ ^ )0= n^ c L x L yq \e^T + 1 (ID 

Using an ansatz of a formal sipjlarity between a pure type II superconductor in high magnetic fields and a current- 
carrying superconductor MakitZI arrived at the following expression for the ratio as/ gin between longitudinal atten- 
uation in the mixed state and in the normal state: 

where A 2 = (|A(r)| 2 ) is the real space average of the gap. As Maki used a semiclassical approximation he did only 
calculate the O'th harmonic. To compare with this result we expand Eq.(0) to first order in A(k). We get: 



^ = l-» (13) 

Oi TV ZKbI 

where (|A(k)|)k is the k-space average of |A(k)|. We see that our theory for the O'th harmonic, which is exact 
within the diagonal approximation produces a term linear in A for the ratio as/a^ as does Maki's conjecture. 
This linear term is somewhat surprising as the Gor'kov expansion of the Green's functions would seem to imply 
that the first correction term is quadratic in A. However, even in the zero field BCS-state one obtainstS as/ot^ = 
2/(exp(A/fcsT) + l) which cannot be expanded in A 2 . So our theory confirms Maki's ansatz of a linear correction term 
for the damping. Thus the non-perturbative linear term in the quasiparticle energy Eq.(|^) coming essentially from 
the degeneracy of the Landau levels shows up in the attenuation result, making a calculation based on the Gor'kov 
expansion questioHable. This is in contrast to the case of the expansion of the thermodynamic potential where contrary 
to earlier claims,t3 it can be proven that these non-perturbative terms cancel, thus making the Gor'kov expansion of 
the thermodynamic potential correct.Ej Our result substitutes Maki's real space average A = ((| A(r) l 2 )) 1 ^ 2 with the 
k-space average of |A(k)|. 

For conventional superconductors in the Meissner state we know that the finite gap for all k supresses the attenuation 
by a factor 2/(exp(A/fcsT + l). We expect that the existence of gapless points will now change this result significantly. 
We therefore split the MBZ into two qualitatively different regions: the "gapped" region where we assume A(k) = 
(A(k)) where (A(k)) is the k-space average of A(k), a the "gapless" region where we assume A(k) = jk v . We 
furthermore assume that this latter dispersion law holds for all A(k) ^ ksT. This is a slight generalization of the 
model used by Dukan and TesanovicQ in their theory for the dHvA oscillations. We assume that the gapped region 
takes up a fraction T of the MBZ. The contribution a(q, w)o, ga p to the attenuation from the gapped region is then 

a(d, u>) Q , gap = —T e(A{k))/kBT + i (14) 

where h = h2ir. As expected the attenuation from the gapped part of the spectrum is strongly supressed due to 
the 2/(exp((A(k))/fcgT) + 1) factor. This result for the O'th harmonic of the attenuation from the gapped part 
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of the spectrum in the mixed state is the same as for the total attenuation in the Meissner state of a conventional 
superconductor. The qualitative new feature comes from the presence of the gapless points. Solving the 2-dimensional 
k-space integral in Eq.(pT|) using the assumed dispersion law, we obtain the contribution a(q, u) g i from the gapless 
part as 

, s ma; 2 fk B T\ 2/n 1-2 1 - 2 /" (2\ (2\ , s 

a(q,« , fl i = Qgi- 2 — r - C - 15 

where Q g i is the number of gapless points in the MBZ. Here C( x ) is Riemann's zeta function and T(x) is the gamma 
function. For a linear dispersion around the gapless points we especially obtain 

mu 2 ir(k B T) 2 

a((l,uj) gl = Q g i——2 (16) 

12ft UJ c qj z 

The relative size of the contributions from the gapped and gapless parts of the spectrum is determined by (A(k)) /k B T, 
T and 7. For (A) lZj"ik B T the contribution from the gapped part can be ignored and the attenuation is given by 
Eq. (|l5|) . Since the normal state attenuation is given by a(q, i^)o,N = m 2 uj 2 /h 2 qh we get 



a(q,u)o,N muj c \ 7 / r\ \rj J \rj 

This calculation is valid in the clean case where the momentum is conserved during the scattering process. As a first 
approximation to the dirty limi|_we can assume that there is no momentum conservation in the scattering process and 
take k z and k' as free variablesEJ (i.e k' z 7^ k z + q). As in the Meissner state, this relaxation does not alter the result 
stated in Eq. ( J17| ) . This is due to the fact we are looking at energies close to the Fermi level such that we can assume 
that the normal state density-of-states is constant. It has been shownEJ that for small impurity concentrations and 
weak scattering potentials the density of states behaves essentially in the same way as for the pure case. We therefore 
believe that our results are somewhat insensitive to the presence of impurities. 

The Landau level structure of the normal state electron energies implies that there will oscillations in the attenuation 
as the external field varies. The oscillatory terms are given by the j 7^ terms in Eq. (pp|) . We end up with the following 
expression for the first harmonic a(q, to) g i.i of the attenuation from the gapless part of the spectrum: 

2ttluo 2 f°° f°° 

a(q,u;)g L1 = Q gi — cos(2n(n f - e(k* z )/hu c )) kdk dx 
nk B l <77r Jo Jo 

x (^ e -V x2h2u i+y 2k2v / k ^T + jj-i^y^fi^J+fi; 2 '! /k B T + j^-i cos(27ra):z;a; e 

/ X 2 h 2 UJ 2 +7 2 k 2 ri 



This integral can be solved in the case of a linear spectrum around the gapless points (i.e 77 = 1). In this case we 
obtain for the 1st harmonic of the attenuation 

k B T x 2 



a(q, uj) g i.i = Q g i —2 cos(27r(n/ - e(k* z )/hui c )) 

qfi u} c n V 7 / 

' tt 2 cosh(2ir 2 k B T/tiw c ) 1 1 

sinh 2 (27r 2 /fc B T/ftw c ) (2Trk B T/hu> c ) 2 \ ' [ ' 

For low temperatures this goes as T 2 . In the case of a general dispersion law around the gapless point given by 
A(k) = jk v , the leading term for the 1st harmonic goes as (T/"/) 2 / 71 . This should be contrasted to the Meissner state, 
or the contribution from the gapped part of the spectrum, where the oscillatory terms will again be damped by a 
factor 2/(exp((A(k))/fcsT) + 1). In the case of the normal state, the 1st harmonic is 

( \ m2uj2 4/cbTtt 2 ru*\/% \\ ron-> 

a q, u) NA = —3 —————— -cos(2TT{n f - e{k z ) hu c )) (20) 

hqti ijj c smh(27r z k B T/?uv c ) 

which is independent of T for low temperatures. In the case of a dispersion law along the z-direction given by e(fc 2 ) = 
t cos(fc z a z ) we have to substitute raj q by 2/ta z \ sin(fc|a z ) — s'm((k~* +q)a z )\ in Ecj.( ^5|) and m cos{2n{nf — e(k*)/fiLo c ))/q 
by 2 cos(27r(n/ — e{k* ^ /Tiuj c )) /ta z \ sm(k*a z ) — sin((fc* + q)a z )\ in Eq.([l^)-(|2C|) where k z ^ are the two solutions 

of the equation t cos((k z + q)a z ) — tcos(k z a z ) = hui. 
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In Fig. 2 we show an example of the acoustic attenuation calculated numerically for two different coupling strengths 
as a function of n/. We have solved the BdG equations self-consistently in 3D as a function of the external magnetic 
field at constant chemical potential. In this paper we have chosen to work with a constant chemical potential. The 
difference between holding the chemical potential constant and holding the number of particles constant is negligible 
in 3D for the normal state.II 2 ] Even in 2D it can be .shown that the superconducting order tends to suppress the 
difference between the two cases in the mixed state.Q The attenuation is calculated via Eq.(||). We have chosen 
parameters such that lud/l^c = 5, kBT/huj c = 0.05 and uj/lo c = 0.01 when nj — 12, where ujd is the usual cutoff of 
the pairing interaction around the Fermi surface. The solid curve is the normal state attenuation which is continued 
into the mixed state to facilitate comparison with the mixed state attenuation; the dashed curve corresponds to the 
coupling strength g/hui c l 3 — 7.85 while the dash-dot curve corresponds to g/h,ui c l 3 — 8.7. The dispersion law along 
the 2-direction is e(k z ) — tcos(k z a z ) where t/hu) c — 0.4. In Fig. 3 we have plotted the corresponding order-parameter 
which in the lowest Landau level approximation for the pairing (which is valid close to H C 2}^ can be characterized 
by the dimensionless number Ao oc A(k). The cowiection between Ao and A(k) or A(r) is given by an obvious 
generalisation to 3D of the results by Norman et aLa The phase transition between the normal and the mixed state 
occurs for u/ ~ 7.6 with g/hui c l 3 = 8.7 and for nf ~ 9.5 with g/hivj 3 = 7.85. As can be seen from Fig. 2, the 
oscillations of the attenuation due to the Landau level quantization persist into the mixed state, although they are 
damped as compared to the normal state oscillations. Eventually the oscillations die out when the off-diagonal pairing 
becomes dominant and the diagonal approximation breakes down. This happens for n/<Lll with g/fiuj c l 3 = 8.7. 
Comparing the attenuation for the two coupling strengths in the region 10 < rif < 11 we get ao(7-85)/ao(8-7) ~ 1.9 
and ai(7.85)/ai(8.7) ~ 2-4 where a(g)t is the i'th harmonic of the attenuation for the coupling strength g. Since 
(A(k)) can be calculated!] directly from A , and the coefficient 7 in the dispersion law around the gapless points 
is proportional to <?A , we can compare the numerical results with the analytical predictions outlined above. If the 
quasiparticle spectrum is essentially gapped, Eq.([l4]) predicts that ao(7-85)/ao(8-7) = /((A7. 85 ))//((A 8 .7)) ~ 8.8. 
If the attenuation is primarily originating from gapless points in the quasiparticle spectrum, Eq.([l5|) predicts that 
a (7.85)/a (8.7) = (8.7A (8.7)/7.85A (7.85)) 2 /''. This gives 2.5 for 77 = 1 and 1.6 for 77 = 2. Hence the numerical 
calculation (ao(7.85)/ao(8.7) ~ 1.9) imply (i) that the gapless points dominate the attenuation in agreement with our 
analytical results and (ii) that the dispersion law is somewhere between 77 = 1 and 77 = 2, since the gapless predictions 
agree reasonably well with the numerical results while the gapped predictions are qualitatively wrong. We cannot 
however make a quantitative numerical determination of the dispersion law around the gapless points. This is due to 
the fact that in order to reduce the computation load, which is high in this 3D case, we have to choose a k-mesh with 
a rather large spacing between the points, (the mesh consists of 100 x 50 points) This means that the gapless regions 
in the quasiparticle spectrum are only probed by a few k-vectors, thereby prohibiting a quantitative determination 
of the dispersion law. Likewise, the analysis above predicts that ai(7.85)/ai(8.7) equals 2.5 and 1.6 for 77 = 1 and 
?7 = 2 respectively. The 77 = 1 prediction of 2.5 agrees well with the numerical result ai(7.85)/ai(8.7) ~ 2.4. Again, 
any quantitative comparison would require a much finer k-mesh. A further complication is that due to the number 
of Landau levels participating in the pairing (approx. 10), the oscillations in the attenuation are quickly damped 
and the diagonal approximation is only valid over relatively few oscillations. This problem could be avoided if we 
were to perform the calculations for experimentally realistic parameters, where there are many more Landau levels 
in the pairing region; however we have not been able to run the programs for such parameters due to the intense 
computation load. The above example does show however that the numerical calculations support our analytic theory. 
In short, the normal state oscillations in the attenuation continue into the mixed state and the damping is dominated 
by gapless points not too far into the mixed state. 

Hence we have calculated the 0'th and the 1st harmonic of the acoustic attenuation in the mixed state. The 
presence of gapless points enhances the acoustic attenuation above the conventional value for the Meissner state. 
When (A) 3ksT such that we can ignore the contribution from the gapped part of the spectrum, the temperature 
dependence of the attenuation is a power law given by a oc T 2 / 1 '. Furthermore, we predict that one should observe 
oscillations in the signal as the external field is varied. The magnitude of these oscillations should have the same 
temperature dependence as the average value of the signal. Hence, by looking at the temperature dependence of the 
attenuation one should be able to detect the presence of the gapless points and to extract the dispersion law around 
these points. From the diagonal approximation it follows that 7 oc (A) where (A) is the k-space (or real space) 
average of the gap. Hence close to the upper critical field H C 2 we expect 7 oc y/l — H/H C 2. This would mean that 
a(q,,u) oc (1 - H/H d2 )- 1/r i close to H c2 . 
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IV. LOW TEMPERATURE 



We will now consider the limit where fcgT <C Tiuj. In this limit only the third delta function in Eq.(^|), which 
describes the creation of two quasiparticles will contribute to the damping. We first calculate the O'th harmonic of 
the attenuation. Using the Poisson identity, making a substitution of variables and using the fact that v s = ui/q -C Vf 
(vf is the Fermi velocity) we obtain the O'th harmonic of the attenuation, when e(fc z ) = k 2 /2m, as 

, s m ^ r~ A(k) r , N X1 e(lu-e) + a 2 

° = ,3 ; ; E / dE[l-f(E)-f(u-E)} [ 1 (21) 

In the limit uj/ksT — > oo this integral can be written as a complete elliptic integralEll and we obtain 



a(q, w)o = <* 2 k£ ^V 1 - 4A(k)7& 2 w 2 J • ( 22 ) 

Here £(k) is the complete elliptic integral of the second kind. The existence of the gapless points again gives rise to a 
qualitatively different result for the ultrasonic attenuation in the mixed state as compared to the Meissner state. This 
is most easily understood by the observation that there will always be attenuation for any frequency in the mixed state 
since there are always quasiparticle states with A(k) < huj/2. Hence the phonon will always have enough energy to 
create two quasiparticles. This is in contrast to the Meissner state, where there is no attenuation for hui < 2A.EJ Thus 
a direct experimental signature of these gapless points would be the absence of the discontinuity in the attenuation 
which is present in the Meissner stateEil, when hu> = 2A, and the presence of acoustic attenuation in the mixed state 
as oj — * 0. We again assume that the dispersion law around the gapless points to leading order is given by A(k) = 7/^' 
(k = |k|) in the region that contributes to the attenuation (i.e for A(k) < hoj/2). Using this we obtain 

a(q '" )0 = Q *W^(2 7 )W XdxEWl - x2V)=n 2^ cq {2^ (23) 
Since the attenuation for the normal state is a(q, w)o,jv = m 2 lo 2 / 271?^ q we have 

-qJ^) VV ^- (24) 



a(q, o;)o,jv \ 2 7/ mu} . 



For the same reasons as for the hui <SC fcsT case, we expect Eq.(|24|) to also be valid in the dirty limit. The remaining 
integral in Eq. ( p3|) can be solved for various 77. We obtain, for instance, I\ — 2/3, I2 — 7r 2 /16, Ii/ 2 = 32/45 etc. 
Again the attenuation has an oscillatory behaviour as a function of the magnetic field H^ 1 . From the Poisson identity 
we get for the 1st harmonic 



a(q,w)i =Q g i^P^ — Y e 2^{n s -e(K)/K^) f kdk f d ^'5(E + E' -u) 
Ah cu c q2TT x J J 

A(fc) 2 



(!-£/£) (!+£'/£') + 



EE' 



e ij6/Kuc t (25) 



We have not been able to solve this integral exactly. However, in the region where the diagonal approximation holds 
(i.e for oj ^ hu) c /A) one can expand the factor e tl i/ huJ <= to a good approximation. This immediately yields that, to 
leading order in uj/ll> c , the amplitude of the first harmonic varies as oj 2 / ri+2 . The next correction term will go as 

So we see that the existence of the gapless points in the MBZ implies that there is a finite attenuation for any 
frequency of the sound wave. There will be no discontinuity in the attenuation as a function of the sound wave 
frequency. As the external field is varied one should observe oscillations in the attenuation. The dependence of 
the attenuation on frequency is algebraic and the power-law is determined by the dispersion law around the gapless 
points. If A(k) = jk 71 we obtain a oc lu 2 /^ 2 . The absence of the discontinuity and the frequency dependence of the 
attenuation should in principle provide the possibility of experimentally determining the existence and dispersion law 
for the gapless points. Again we expect a(q, u) oc (1 — H / H^)^ 1 ^ close to H C 2. By making the same substitutions 
as in the Kuj <C fc^T limit, one can obtain the results for the case e(k z ) = tcos(k z a z ) relevant for layered structures. 
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V. CONCLUSION 



In this paper we have considered the acoustic attenuation in a type II BCS superconductor at high magnetic fields 
using both numerical and analytical methods. We have shown that away from the semiclassical regime where the 
Landau level structure of the electronic states is important, the attenuation will in general be an oscillatory function 
of the external magnetic field. Furthermore, since the attenuation probes the quasiparticle density of states, the 
presence of gapless points in the quasiparticle spectrum makes the attenuation qualitatively different as compared to 
the Meissner state attenuation. For fc^T -C hu) c the attenuation is an algebraic function of the frequency and there is 
no discontinuity as opposed to the Meissner state attenuation. For fc^T 3> Tiuj c the attenuation is an algebraic function 
of the temperature. The exponent of the power law is determined by the dispersion law around the gapless points. 
This behaviour should in principle be experimentally detectable; such an experiment would provide confirmation of 
the existence and nature of the gapless points. 
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Figure Captions 



Fig. 1: The lowest quasiparticle energy in units of Uto c for two different values of the order parameter ((A(k)) ~ 0.3hui c 
and (A(k)) ~ 0.057iw c ). The k-vector is measured in units of 2n/L x where L x is the size of the sample in the 
cc-direction. The solid lines show the exact numerical result while the dashed lines correspond to the diagonal 
approximation. 

Fig. 2: The attenuation as a function of rif. The solid line is the normal state attenuation. The dashed line is for the 
coupling g/huj c l 3 = 7.85. The dash-dot line is for g/hu)J 3 = 8.7. 

Fig. 3: The order-parameter as a function of n/. The dashed line is for the coupling g/huij 3 = 7.85 while the dash-dot 
line is for g/hui c l 3 = 8.7. 
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Figure 1 (Bruun, prb) 
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